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FOREWORD 


This document has been prepared jointly by the Facilities Support . 
Office of the Earth Observations Division , Lyndon' B. Johnson Space 
Center (JSC), Houston , Texas,, and by the Computer Science De- 
partment of the Space Information Systems Operation ( SISO ), Aero- 
nutronic Ford Corporation. It was published by SISO in accordance 
with EO-004N, Crop Assessment Subsystem (CAS) data system, of 
Contract NAS- 9-1261, Schedule V. 
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SECTION 1 
INTRODUCTION 


1.1 PURPOSE 

This document describes the functional data processing requirement 
for the Crop Assessment Subsystem (CAS) of the Large Area Crop 
Inventory Experiment (I.ACIE). These requirements shall be used 
as a guide for software development and implementation. The soft- 
ware implementation plan will be described in a detailed design 
document to be published at a later date. 

1 . 2 SCOPE 

This document is limited in scope to requirements pertaining to 
the production processing software of the CAS. In general, these 
requirements include the following: 

• Command input processing 

• Processing and storage of data from external sources 

• Wheat area, yield, and production estimation 

• Computation of standard statistics associated with the crop 
estimates to assess the reliability of the estimates 

• Aggregation of crop estimates to the various hierarchical 
elements of the active LACIE countries 

• Generation of aggregation reports. 

Paragraph 1.4 of this document presents a brief overview or tne 
CAS software system components. 

Section 2 describes the CAS/anaiyst interactions and defines the 
support software routines. Section 3 details the CAS applications 
software; Section 4 provides CAS data base requirements; and Sec- 
tion 5 describes Classification and Mensuration Subsystem (CAMS) 
process ing . 
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1^3 APPLICABLE DOCUMENTS 

The following document has been used for support data in the 
preparation of these requirements. It may be referenced by the 
reader for obtaining additional information in specific areas 
of .interest. 

Large Area Crop ' Inventory Experiment (LACIE) Level III Base- 
line Crop Assessment Subsystem (CAS) Requirements 3 LACIE- 
00200, Vol . IV, 16 December 1974.. 
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1.4 SYSTEM OVERVIEW 

The CAS software system shall be implemented in the multi-user 
environment provided by the Digital Equipment Corporation (DEC) 
RSX-11D Operating System on the PDP 11/45 Computer System in the 
DTL. The software may be exercised in both a batch mode and an 
interactive mode of operation. Each user of the CAS software will 
operate from either the master data files for a country or a work- 
ing copy of the data files for that country on a non-interference 
basis. Access to the CAS software and its data files shall be 
restricted to authorized personnel via the standard RSX-il!) user 
identification code (UTC) and password security features. 

The CAS master data base is made up of multiple data files, one 
for each type of data for each active LAC1E country for each crop 
year. The four data types are: 1) allocation data, 2) classifi- 

cation results, 3) wheat yield data, and 4) historical agricul- 
tural data. Each data type? corresponds to an external data source. 
The allocation data shall be utilized to define the logical struc- 
ture of the remaining data files and may not be updated on a cas- 
ual basis. The remaining data files shall he routinely updated 
by data base in i t i a 1 i z at i on/preproces s ing software as new data 
becomes available. Only the classification results data shall be 
preprocessed to yield an intermediate product. 

Each application program within CAS is designed to operate on the 
data files for a single country. Thus, each CAS user must estab- 
lish his definition of data files with respect to a single country 
and whether the master data files or a working copy of the data 
files arc to he utilized. The data base management software shall 
then provide the means for the applications programs to retrieve, 
update, and store data. 

The user may select tc execute one of the following processors: 

• Area/yield/production est inat ion/aggregat ion 

• Report generator 

• Data base change programs. 
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Control inputs shall be accepted, validated, and processed by 
either the interactive input command processor or by the batch 
input processor. After control is passed to the selected proc- 
essor and processing is completed, output shall be generated 
via the display processor, a subset of the report generator. The 
relationship between the various software components of the CAS 
is illustrated in figure 1-1. 


1-4 



JSC- 10009 



Figure 1-1 CAS Functional Data Flow 
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1.5 REVISIONS 

Because this is a controlled requirements document, any required 
changes shall be submitted to the chairman of the Facilities 
Change- Control Panel for review and approval. 



CAS SYSTEMS 
SOFTWARE 
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SECTION 2 

CAS SYSTEM SOFTWARE 


2.1 GENERAL 

The CAS system software shall comprise three major areas: inter- 

active terminal software, system utilities, and support software. 
The interactive terminal software shall provide the interface 
between the analyst and the application program. The system 
utilities shall be provided to assist in the implementation of 
the CAS software system. The support software shall be provided 
as a library of routines in support of the application program. 

2.2 INTERACTIVE TERMINAL SOFTWARE 

The CAS interactive terminal software shall provide for terninal 
initialization, input query operations with the operator, and CAS 
report generation and display. Interactive software shall also 
provide operator interface to the Data Base Change Program. 

2.2.1 CAS Input Processor . The CAS input processor shall ini- 
tialize the terminal cathode ray tube (CRT) and process analyst 
input commands. The CAS input processor thus provides the inter- 
face between the analyst and CAS application programs. 

2.2.1. 1 I nitialization . CAS terminal initialization shall be 
performed by the RSX Operating System. To initialize a CAS ter- 
minal, the analyst shall strike the control C key on the CRT 
terminal. This will activate the Monitor Control Routine (MCR) 
which shall display a prompter. The analyst shall now "sign-on” 
to the system by entering the HELLO command and a CAS '.IIC number. 
CAS files shall be protected and may be accessed only by use of 
the CAS UIC. The system shall now request that a password be 
entered, but the password shall net be echoed back to the CRT 
screen. The system shall validate the password. If the code is 
illegal, an error message shall be displayed and another request 
made for the password. If the password is accepted, the analyst 
shall be logged on to the system and may run any interactive CAS 
programs as well as other system programs. 
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2.2. 1.2 Input Command P r oces sor. To initialize or activate the 
CAS software system, tile CAS analyst, shall type CAS after the 
MCR prompter. This shall cause the CAS input command processor 
to be loaded into memory ' The command processor is new ready 

to accept input commands from the analyst. A command menu (list) 
shall he displayed for the analyst to select the desired function 
by typing the appropriate three letter command name. 

After the command is selected, the input command processor shall 
begin to request processing parameters by question and answer 
interaction with the analyst. The input command processor shall 
output a message to the analyst requesting a particular input 
parameter. A list of the appropriate input options may also be 
displayed for the analyst. The analyst shall transmit the input 
parameter back to the computer by typing the parameter on the 
next line. The input command processor shall validate the input 
and either display an error message or continue to the next param- 
eter request. After all parameters are collected and validated, 
a run request shall be made to the selected application program. 
The input processing parameters shall be transmitted to the called 
program. Appendix A discusses a typical operator/CRT interaction. 

Upon completion of the application function, control shall he 
passed to the display processor (sec paragraph 2.2.7) to output 
the results of the application program. The input command proc 
essor shall or.ee again receive control upon termination of dis- 
play processor output. At this time the analyst may request 
additional processing or exit back to the operating system. 

2.2. 1.3 Termi na t i o n . After the display processor has completed 
output of all requested reports, control is returned to the input 
command processor where the function or command menu is displayed. 
One of the commands available to the operator shall be an exit 
command; e.g., to terminate the CAS system, the operator shall 
type the three letter command TXT for exit. This shall return 
control to the RSX Operating System. At this point, the operator 
signs -off the terminal by typing the command BYM which completes 
the termination procedure. The only valid input command at this 
time is the HELLO command which is the first step of initializa- 
tion (see paragraph 2. 2. 1.1). 
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2.2.2 - Display Processor . The display processor shall standardize 
the techniques for generating a CAS report. The display processor 
is. a user-oriented service program which shall generate a display 
based on input parameters supplied by the user. Multiple output 
devices shall be supported by the display processor; i.e., the 
reports may be displayed not only on the analyst's CRT terminal, 
but also may be output to the line printer and/or magnetic tape. 

The output devices to be used may be in any combination and shall 
be required input parameters to the display processor. 

• 2.-2. 2.1 Modes of Operation . The display processor shall operate 
in either automatic or manual mode. In the automatic mode, only 
the line printer and/or magnetic tapes shall be utilized. Once 
started, the display processor shall output the requested report (s) 
to the appropriate device (s] without analyst intervention. 

The manual mode is implied when the CRT is requested as an output 
device for the CAS reports. In the manual mode, the summary report 
shall be output to the CRT terminal and any other devices requested. 
The display processor then waits for the analyst to enter a com- 
mand. At this point the analyst shall be able to request the next 
sequential report, to request a different report, or to exit the 
•display processor (see paragraph 2. 2. 2. 3). An exit request shall 
return control to the input command processor, which once again 
displays the input command menu. 

2. 2. 2. 2 Input Parameters . The CAS application programs shall 
exit, after storing all generated report data on disk, by schedul- 
ing- the display processor. Two input parameters may be transmitted- 
.to the display processor at this time. The first parameter shall 
be the output device or devices. Valid output devices for CAS 
reports sha.ll be, the CRT terminal, line printer, magnetic tape, 

or any combination of these. The second parameter shall be an 
additional reports indicator and shall only be valid with the line 
printer and magnetic tape devices. This parameter shall also in- 
dicate that all individual reports are to be output in addition 
to the summary report which is automatic. This- parameter shall 
not apply if the CRT output device has been selected, as the oper- 
ator has the opportunity to request individual reports after the 
summary report has been displayed. 
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2. 2. 2. 3 Display Commands . If the analyst selected the CRT ter- 
minal as an output device for the CAS reports, the display proc- 
essor shall wait for an analyst response after output of the sum- 
mary report. The analyst may respond with any valid display 
processor coimnand. The three basic command functions are dis- 
cussed in the following paragraphs. 

A. Report Retrieval . To retrieve a CAS report, the analyst 
shall input one of the report retrieval commands. Each 
command shall be composed of a 3-letter command code and 
a report key number. The first letter (A, Y, or P) of 
the command code shall indicate whether the requested 
report is an area, yield, or production report. The. next 
two letters shall indicate the report level to be re- 
trieved. The following options shall be available: 

• Summary 

• Region 

• Zone 

• Strata 

• Substrata. 

The report key number is a 4-digit integer indicating the 
hierarchy element number. Each report contains the key 
numbers of hierarchy elements of the immediate lower level 
which were used to produce the report. (For explanation 
of CAS .reports and report contents, see paragraph 3.5.) 

B. Exit . To exit from the display processor, the analyst 
shall input the exit command (EXT). This command does 
not require any additional input parameters. The exit 
causes return to the CAS input command software. 

C. Device Change . The device change command (DEV) shall allow 
the analyst to select a new device or combination of devices 
to be used for output of CAS reports. The analyst inputs 
the command DEV along with the device mnemonics. All sub- 
sequent requested reports shall be output to the new de- 
vice (s) . 
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2.2.3 Data Base Change Interactive Software . The data base 
change interactive software shall provide the interface between 
the operator and the Data Base Change Program. The interactive 
software shall accept input commands and display the results of 
the processed command. ' ' 

2. 2. 3.1 Initialization . The data base change interactive 
software shall initialize the CRT terminal by displaying a re- 
quest to enter a command. Retrieval and exit shall be the only 
valid commands at this time. 

2. 2. 3. 2 Data Base Change Command Processor . The data base change 
command processor shall accept command inputs from the analyst. 

The commands and any parameters shall be validated and the in- 
formation passed to the Data Base Change Program for processing. 

If an invalid input is detected,' an appropriate error message 
shall be output along with a request' to re-enter the command. 

An input command prompter shall be used to request an input from 
the analyst. The prompter shall display the valid input commands 
and the command codes. Command codes (see paragraph 2. 2. 3. 3) 
shall be unique alpha characters. 

2. 2.3.3 Input Commands . The command processor shall accept five 
commands from the analyst. The following is a discussion of each 
of these commands. 

A. Retrieval (RET) . The retrieval command shall be used to 
retrieve requested data base information. This information 
shall be formatted and displayed on the analyst CRT ter- 
minal. Retrieval of the data base information is a pre- 
requisite to performing a data base change. The re- 
trieval command shall be a valid input at any time. Fol- 
lowing-input of the RET command, the command processor 
shall request input of parameters specifying the exact 
data base information being requested (see paragraph 4.4). 

B. Change (CHA) . The change command shall be used to input - 
data base change information. The CHA command shall be 

a valid input following the RET command. The change in- 
formation shall be displayed along with the information it 
is to replace. CHA does not cause the CAS data base to be 
updated; the actual change of the data base is accomplished 
by the store command. 
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C. Store ( STR) , The store command shall be used to write the 
change information to the CAS data base. The STR command 
initiates the update causing all change information to he 
written to the CAS data base. The command may be entered 
following a CMA command. After the command is entered, 
the command processor shall request store information from 
the analyst. This information shall include the reason for 
the change ar.u the name of the analyst. 

D. Exit (HXT ). The exit command shall be used to exit the 
Data base Change Program. Control shall he returned to 
the CAS input command software for subsequent display of 
the CAS command menu. 

K. Page (PAG) . The page command shall he used to request the 
next page of a multipage display. This command shall be 
valid following the RhT command. No additional parameters 
shall be required. Additional inputs of the PAG command 
shall cause subsequent pages of the display tc be output. 

F . B ack P age (PPG). The back page command shall allow the 

analyst to request the previous page of multipage display. 
This command shall he valid following the PAG command, and 
no additional parameters shall he required. 
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2.3 SYSTEM UTILITIES 

The following system utilities shall be converted to augment those 
utilities provided by the RSX-11D Operating System. 

* v ' 

2.3.1 Nonstandard Magnetic Tape Copy . This program shall copy 
an 800/1600 bpi, 9-track magnetic tape from one tape drive to 
another. The user shall specify the input and output tape drives. 
Optional parameters may be input to position the tape by file and 
record. The processing shall start by reading a block and check- 
ing for errors. If an error is detected, an appropriate message 
shall be printed with the block number and block size. Whether 
an error was detected or not, the record shall be written to the 
output magnetic tape. If an end of file (EOF] is detected, a 
message shall be written to the line printer, and an EOF written 
on the output tape. A check shall be made for three consecutive 
EOF's. When detected, a message shall be written to the line 
printer and processing terminated. 

2.3.2 Nonstandard Magnetic Tape Dump . This program shall dump a 
800/1600 bpi, 9- track magnetic tape to the line printer in octal, 
hexadecimal, and ASCII. Parameters may be input to position the 
tape by file and record. The processing for a block shall start 
by reading a record and checking for errors. If an error is 
detected, an appropriate message shall be printed. In any case, 
processing continues by printing the block number and size. This 
shall be followed by a check for EOF. Unless the block is an EOF 
mark, the ‘contents of the record shall be dumped on the line 
printer. If the record is an EOF mark, an EOF message shall be 
printed and a check made for three consecutive EOF's indicating 
logical end of tape (LEOT) . IF LEOT is not detected, the next 
block shall be processed. If LEOT is detected, a message shall 
be printed and processing terminated. 

2.3.3 Tape Merge . This program shall merge a single file of a 
800/1600 bpi, 9- track tape to a common merge tape. The merge 
program shall write the single file at the end of the merge tape 
and write three EOF's after the file. Error conditions shall be 
checked after each read and write. If an error occurs, an error 
message shall be printed and three EOF’s written on the output 
tape and processing terminated. A summary message indicating 
total number of files on the merge tape shall be output at termi- 
nation. 
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2.4 SUPPORT SOFTWARE 

The following support software routines snail he provided. 

2.4.1 Byt e Sw ap pin g ( BSWAP) ♦ BSWAP shall be provided to allow 
swapping the position of the bytes in words in an array. The ad- 
dress of the array in which the bytes will be swapped and the word 
length of the array shall he passed to BSWAP in the calling se- 
quence. BSWAP shall swap the bytes in each word in the array, 
beginning with the first word in the input. The byte-swapped 
words shall be stored in an output array whose address is passed 
to BSWAP. 


Byt e Ar ray Move f rom P acked Arra_v ( BMOV) . BMOV is 
mtir.e which shall be provided to unpack 


sembly language 


an a s - 
d a t a 

from one array (two bytes per word) into another array (or.e byte 
per word). The destination array must be dimensioned twice the 
source array. The high order byte of each word in the destina- 
tion array shall be cleared. The number of bytes to be unpacked 
shall he an even, number . 


2.4.3 Bvte Array Move, from. ; J n p a c k e d Arrav (WVOV). WMOV is an 
assembly language routine which shall be provided to pack data 
from, one array (one byte per word) into another array (two bytes 
per word). The number of bytes to he packed shall be an even 
number. The dimension of the source array must be twice that of 
the destination array. 


2.4.4 Array Move (AMOY) . AMOY' is an 
which shall be provided to facilitate 
data very rapidly. The input paranet 
source array, the destination array, 
be moved . 


assembly language routine 
the moving of blocks o: 
ers to AMOY shall be the 
a n d t h e number of w o r d s to 


2.4.9 
wh i ch 
accept 
the N 
a real 


Zeros Array (Z ERO) . i 
shall be provided to ze 
a starting address and 
16 -bit words beginning 
a r r a y, N must he twice 


FRO is an assembly language routine 
ro an array. The routine shall 
a count N. ZERO shall sec to zero 
a t t h e s t a r f i n g add r c s s . To clear 
the number of real values to be 


c 1 eared 
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2.4.6 Bit Extraction (I F LO). IFLI) shall extract the desired 
bit string given the starting hit number, the number of bits to 
extract, and the 16-bit word to be accessed. The bits shall he 
numbered 0 through 15, counting left to right. The function 
result shall be returned ‘right -justified , and remaining bits set 
to 0. 

2.4.7 Bit Inser t ion (1SLT). I SET shall be used to set bit 

strings into a lo-bit word. The user shall be required to pass 
the starting bit number and the number of bits to set in a given 
word. 1SKT shall be a function subroutine and shall not alter the 
contents of the calling arguments. 

2.4.8 J ulian Calendar Hate C onversio n CniATL) . JDATE shall pro- 
vide data conversion which gives the user three options of con- 
version. These options are: 

• Input a Julian date and output the Gregorian date and the 
Julian day of the year 

• Input a Gregorian date and output the Julian date and the 
Julian day of the year 

• Input the Julian day of the year and the Gregorian year and 
output the Julian date and the Gregorian month and day. 

2.4.9 Formatted Tape Dunn . Tape dumps shall be provided to 
dump the following tapes: 

o Classification and Mensuration Subsystem (CAMS) /CAS inter- 
face tape 

• CAMS input tape 

• Yield Estimation Subsystem (YKS) input tape 

• Data Acquisition, Preprocessing, and Transmission Subsystem 
(i)APTS) input tape. 

The user shall have the capability to select specific records and/ 
or files on the input rape. 
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2.4.10 Ma trix M u ltip lic atio n (M AT.V'J l.) . MATMtJI. shall provide 

the capability to multiply an array whose dimension as defined 
on input times another variably dimensioned matrix or vector. 
The routine shall compute any one of the fol lowing, depending on 
an option indicator: ' ' 

• MATRIX * VECTOR 

• MATRIX TRANSPOSE * VECTOR 

• MATRIX * MATRIX 

• MATRIX TRANSPOSE * MATRIX 

• MATRIX * MATRIX TRANSPOSE 

• MATRIX TRANSPOSE * MATRIX TRANSPOSE. 


Inputs to MATMtlL shall be: 

• Input MATRICES (2) or MATRIX and VECTOR 

• Output MATRIX or vector location 

• Option code indicator. 

2.1.11 Statisti cal Error E ur.c ii on (ERE). ERE is a routine which 
shall be provided to approximate the statistical error function 
that is mathematically defined as: 

7 

? /*X . T ~ 

e r l i x ) -■ / e d z 

v- Jo 

2.4.12 M a t. r i x ■ r. y c r s i o n ( M A T I N\‘ ) . MATINV is a routine which 

shall be provided to determine the inverse of an N> N matrix M 
for N < !>(). The inverse matrix, M l , is defined such that 

M M' 1 ■■ I 


2 - 10 



JSC-10009 


where I = the identity, matrix 


10 0 
0 10 
0 0 1 - 


0 

0 

0 


0 0 0 


1 


2.4.13 Conversions (A2EBCD/EBCD2A) . A2EBCD and EBCD2A are con- 
version routines which shall provide conversion from ASCII code 
to EBCDIC code and vice versa. The routine's shall have the follow- 
ing calling- parameters : 

• Input array to be converted 

• ‘ Output array that has been converted 

e Number of characters to be converted (INPUT and OUTPUT) 


« Status of the conversion (error detection) 


2.4.14 Metric/USA Conversion . Metric/USA units conversion shall 
be accomplished by multiplication with standard conversion con- 
stants. The following is a list of the conversion constants. 


1 Acre 
1* Hectare 
1 Bushel (wheat) 

1 Metric Ton 
1 Bushel/Acre 
1 Quintal/Hectare 


= 0.40468564 Hectares 
= 2.4710538 Acres 
= 0.02721552 Metric Tons 
= 36.743740 Bushels 
= 0.67251015 Quintals/Hectare 
= 1.4869664 Bushels/Acre 
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SUCTION 3 

CAS APPLICATIONS SOFTWARE 


3.1 GENERAL 

The CAS applications software shall provide the capability for 
area, yield, and production estimates with standard statistics 
for spring and winter wheat, for each stratum, r.onc, region, and 
country in the LACIK. Also, the applications software shall be 
capable of generating reports of area, yield, and production 
estimates and statistics. Units for all computed parameters 
shall be in the metric system of mcasureir.ent. 

The software shall support batch or interactive operations by the 
CAS analyst. The analyst shall have the capability to permanently 
redefine the correlation between the elements of the sampling 
strategy hierarchy for region membership only. 

3.2 AREA ESTIMATION AND AGGREGATION 

The area applications software shall calculate area estimates and 
standard statistics from the input data (paragraph 3.2.4). The 
area estimates and variances shall be aggregated to the level 
specified by the process specification parameters (paragraph 
3. 2. 3. 3). The area and statistical information shall be output 
(paragraph 3.2.5) as requested in standard report formats. 

3.2.1 Area Iis_Mma_r cs . Wheat area estimation shall be performed 
at each defined area level of each country. Separate estimates 
shall be performed for spring and winter wheat. 

The lowest subdivision (base level) of a country shall be either 
strata or substrata in all cases. Each of the elements cf the 
base level shall be classified into one of three possible groups: 

• Group I - Stratified 

• Group II - Probnbi i i t.y- Proper t iona: - to - Si zc (PPS) 

• Group III - Ratio. 

The PPS group shall exist only if the base level is the substrata 
level . 
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3. 2. 1.1 Strat i ( ied Subst ratn/St rata . Stratified area estimation 
shall be performed on a base level of substrata (or strata) ele- 
ments. The wheat area estimate siiail be calculated by the follow- 
ing equations : 

Mjk 

'V - <wv £ A jki (i > 

i»l 


A 


U 


L 


k= 1 



( 2 ) 


Where : 


A 




LI 


A 

N 


ljk 

j k 


Combined estimate of all Croup I subelements (substrata 
or strata) in t b. c j t h strata 

Number of Group I subelements for which classification 
data is available 

Kstimate of the kf-Ji Group I subelenent (substrata or 
strata) in the j t h strata 

Number of agricultural segments (after exclusion) in 
the ktb subelement in the y.h strata 



Number of sample segments for which classification data 
is available in the k th subelement in the j?./: strata 


R 


j* 


Ratio of the true subelenent area to the gross pseudo 
subelenent (before exclusion) area 


A . , . 

JKt 


estimated wheat area in hectares for the it't sample 
segment in the kr/; subelement in the j th strata. 


Any subelener.ts for which 
shall be accounted for by 


no classification data is 
the Group III estimator. 


avai lab 1c 
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3 . 2 . 1 . 2 Probabi 1 i t.y- Proport io nal-to-Sir. c (PPS) Su b s t r a t a . P P S 

area estimates shall be computed tor the group of elements within 
a strata that is designated as PPS elements. The wheat area 
estimate shall be calculated using the following equation: 


A 


2j 



<W 




(3) 


Where : 


A 


jk 


A j y | since only one sample segment may be allocated to 
any PPS subelement 



Kpoch year wheat, area for the kih subclement in the PPS 
group of the j th stratum 


W . 
J 



K . . 
J* 


M 


2 i 


Number of sample segments in the PI’S group of elements 
for which classification data are available 


1.2 *= Number of PPS substrata in j ih strata 


N 


jk’ 



and R.. 

J * 


are as previously defined. 


3. 2. 1.3 T otal Str ata _Area lis t in are . The wheat area estimate of 
each strata shall be computed ns the sum of the Group I, Group II, 
and Group III component estimates which comprise the strata. 

The Group II! component shall not be individually computed, but 
shall be accounted for in the ratio estimation for the strata as 
presented ir. the paragraph 3.2. 1.4. 


A. - A. . + A . ♦ A. . (4) 

J 1 J 2 j o j 


If the strata dees not contain a Group III element, the sum of 
the Croup I and Group II components shall define the strata 
estimate. 


A. 

J 




(S) 
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3. 2. 1.4 R atio Suhs t rata /Strat a . The area estimate for each 
stratum containing a Group III element shall be selected from 
one or more candidate estimates according to a minimum variance 
logic. Each of the candidate estimates shall be computed from 
a common algorithm which employs different input data sets. An 
estimate shall be computed for each hierarchical level element 
that contains t ho Group ill element and also contains one or 
more observations. 


Thus if the jtfc stratum contains observations and a Group II! 
substrata, belongs to the Lth zone, nu/; region, and nth country, 
then four variance estimates shall be computed vising the four 
hierarchical spaces that contain the jth strata. 


V(A. n ) « fjA^, V (A , n ) , 




> ' 



J » 




] 


V (A j^) , V (A : K ) , and V(Ai J ) are similarly defined over the m, 
i , and j spaces, respectively 


Where : 


V(A. ) r Variance of the area estimate for the jth strata 
- computed from data belonging to the space n 

A n n Area estimate o: a (-roup I element belonging to space 


V (A ) - Variance of Aj 


n 


Historical wheat area data of A 


n 


. . n 


o 


A,', V(A, *), and W,' are similarly defined 

- Historical wheat area data for the Group ill component 
of the jfJ; strata. 
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Variance estimates at the zone, region, .and country levels nay 
also be computed from spaces of which they arc elements. 

V(A t n )i V(A £ m ); V(A. £ ) - Zone level 

V(A m n ); V(A m ) ~ Region level 

VCA^ 11 ) ~ Country level 

These variances are not an output but shall be used to select the 
strata estimate space. 

The strata estimxite space (s) shall correspond to either the strata 
( j ) , the zone (C), the region (m) , or the country (n) as determined 


by the following 

logic : 


If V(A 4 l ) 

> 

r.i 

J 

then s — j 

If V (A m ) 
rn 

> 

it n 

v(a/). 

then s - l 

If V(A n n ) 

> 

y. 

men 

V (A "*’) , 
m 

then s ~ n 

Otherw i se 

s 

- n 



Once the 

strata < 

est ir.are 

space (s) has beer, chosen, both the area 


estimate and the variance estimate can be established. 


A. «= A . s 
J J 

V(A.) - V(A S ) 

Exact equations for the strata estimates are: 


A. - A. . ♦ A, . ♦ r . S (A. S +A, S ) 
J lj 2;, j 1 2 ’ 


( 6 ) 
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Where : 

Tj S = EW 3 j /(EW 1 S +2W 2 s ) 

9 

W = Historical wheat area of a Group I element belonging 
to s for the epoch year data 

s 

W 2 = Historical wheat area of a Group II element belonging 
to s for the epoch year data 

Wj^ = Historical wheat area of a Group III element belonging 
to strata j for .the epoch year data. 

3. 2. 1.5 Selection Criteria for Classification Results . Classi- 
fication results shall be rated on a scale of 0 to 15, with the 
higher numbers representing the best results. A threshold number 
can be specified, and the most recent classification results 
which are equal to or greater than this threshold value shall be 
used for area estimations. Other threshold criteria, as spec- 
ified in paragraph 3. 2.3. 3, shall be evaluated in the selection 
process. 

3. 2. 1.6 Dynamic Regrouping of Group III Substrata/Strata . Each 
stratified base level element and each PPS-element group shall be 
dynamically and temporarily redesignated as Group III elements if 
no classification data is available for that element. However, 
if at least one sample segment classification is available, then 
the area for that element shall be computed according to the 
nominal grouping of the element. 

3.2.2 Standard Statistics . The standard statistics for area 
estimation and aggregation are composed of standard deviation (a), 
coefficient of variation (CV) , 90 percent confidence limits (CL ) f 
and probability of achieving 10 percent error or less (fl) . Each 
statistic is a function of variance (V) as follows: 

° 5 - 

Where s 55 Any area element 
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CV » o /A 
s s s 


Where A - Area estimate of the s th element 
s 

• K 
(♦) 


( 8 ) 


Cl. 


A ♦ (1.645) o 
s s 


CL - A - (1.645) 0 

S 5 S 


« s - erf J(l/>/?)!l/(10CV s )]j 


(9) 


( 10 ) 


Where : 

erf (0A ft) 


n Hrror function associated with normal probability 
curve 

- Standardized random variable 


S.2.2.1 St rati fied Area Ksti matc V ari ance . The variance of 
wheat area for stratified (Group I) estimates is given by: 

V "<V " v u ' A V < K jk Z/M jk5 V 01 


Where: 


1.1 " Number of Group I elements for which classification data 

is available i r. the j th strata 

f. L - M . . / N . , 
j k j k j k 


V 


Ratio of the true element area to the gross pseudo cle- 
ment (before exclusion) area for the element of the 
jt/: strata 


N' u Number of agricultural segments (after exclusion) in the 


jk 


M . . 
Jk 


kth element of the )th strata 

Number cf sample segments for which classification data 
is available in the kr/’i element of the jth strata 


i - / 
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S ‘ 
D k 


M 

jk 2 

1/M jk E °Jki 

i=l 

or 

M.. 
jk 


if SC^ = 0 


a >2 * r n /~i ^ 


^V 1} E CA jk rV if SC 1 ' 0 

i=l 


( 12 ) 


If S 2 > C 2 /4, then S., 2 = C 2 /4. 
jk 1 jk 1 


= Number of hectares per sample segment 

A.^ = Estimated wheat area (in hectares) for the itTz sample 

segment, kth element of the j th strata 

A., = Mean wheat area for the sample segments in the k th 

■* element of the j th strata 

SC^ n = System constant for country n 

= 0 for countries with detailed historical data (e.g., 
United States, Canada, Australia) 

f 0 for other countries. 


2 2 -1 

a jki = g jk ” D (the bar below a 

vector) 

[C] = An H x H matrix (SC 2 < H < 50) 

T t -1 2 

c.. = [i - x. 1 (x x) x. ] 

ii 1 i v i J 


C ih = [ X i T (X T X) _1 X h ] 2 ; i = 
— — — h = 

i t 


symbol implies a 

whose elements are: 

, i 1 , 2 , . . . ,H 

1 , 2 , . . . ,H 
1 , 2 , . . . ,H 
h 
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SC 2 = System constant designating the minimum num- 
ber of sample segments over which to perform 
a regression. The nominal value for SC ? will 
be 10. 


H 



+ 


LI 

E V 


k=l 


If H < SC 2 , the computations are done at the 
next highest level where the number of sample 
segments with classification data > SC 2 . 

M- . = Number of sample segments in the PPS 
J group of elements for which classifica- 
tion data are available. 

H is the total number of sample segments 
(stratified and PPS for which classi- 
fication data is available in the small- 
est space containing the H observations). 


X = An H x 2 matrix 


* 

1, x x 


1 

r*H 

1 

1, x 2 * 

= 

T 

V 

« 


• 

• 


• 

• 

* 

_ 1, X H _ 


• 

H J 


x ki* = M N A) for 1 = 

W k = E P ocl1 year wheat area for the k th 
substrata in hectares 
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D = [d^, d 2 2 d H 2 ] 

d = A - x£ 

6 = Cx t x)" 1 x T a 

A ss A.. . 
jki 


3. 2. 2. 2 PPS Area Estimate Variance. The variance of wheat area 

(1*3) 


for PPS (Group II) estimates is given by: 


Var (A.,.) = V 


Where : 

L2 

V_. 

2 l 

■E R 


k=l 

- 

L2-1 

V " 
2 j 

■ E 


(i'EY) 


/n ji ] 


i=l i'*=i+l 

V 2 j > “ 0 if only 1 year of historical data is available 

p = [w (EY) (EY) • 
jk L 3 k /W j J 

H =M P..; P.. = [W./ EY Vw. (EY) ] 

1 1 2j 3 i 31 L 3 1 2 

n . M p -p - rw ( EY )/ W (EY)! 

ji' M 2 j j i" * P ji" [W ji- /W j ] 

EY - epoch year historical data 

^EY ~ the most recent year that is not the epoch year 


(14) 
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K CUV) . £ (E 

J tl Jk 


0:Y) 


n . . . 


= I(M 2j -i)/M 2j ] :: ;i n 3 ., ♦ ( C» 2j - - 1 )/ (* 2j ) 1 


0' ■’ 
U ‘ji ”ji' 'ji'ji 


[(M 2j -D/(M, ; 


7 I- 2 - 

V M . “ 
2j J 1 j 1 fr'i -* 3 


♦ [ 2 (M, . - 1 ) / (m ) ^) cn .. . _*n . .n . 2 n . . }) 

2j 2.1 J1 ji Jiji J1 J1 


(J(H -U/CM ) 4 1 (3 2 n *3 Bji J ) 


|2ft V 1,/( V" ! £n Ja s 


0.-1 




0 V 

Jk 


v _ 


liquation 12 for PS’S substrata with a segment or 
( v i i + Y ■> \ * i v ' 


lj '2; jk' , 


for PPS substrata without a segment 


2 2 ? 2 

If > Cj /4 , then r c ~| /4 , where C f 3 number o 

hectares per sample segment. 
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2j 


H 

£ ‘W - 

h=l 

H 

am £ z jh 

h=l 


aM* 

i-Z i 

£ a. h ) z - ci/hd 

h=l 

' H 

2 X jh 
_h=l 

2 


CSC 2 < H < 50) 


'U ' C1/H) ti Z J h ' tY2 5 /H ’ S V 


CSC 2 < H £ 50) 


Z jh 

V 


*°*e Ca jh ^ 

*°*e p jh + *°*e C1 * p jh ) 


P,- k - [»4- 


CEY) 


'jh ” L "jh for each sa mple segment h in the 

j th strata 


= X_. ^ for PPS subelements. 


If H < SC 2 , the computations are done at the next highest level 
where the number of sample segments with classification data 
H > SC 2 . 


3 * 2 *2.3 Total Strata Area Estimate Variance . The variance of 
the wheat area estimate of each strata that does not contain a 
Group III element shall be the sum of the Group I and Group II 
variance components. 


V. 

3 




(15) 


If the strata contains a Group III element that is independent 
of the Group I and Group II components, the strata variance 
estimate shall be defined by the sum of the Group I, Group II, 
and Group III components. 
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V. 

J 


V 


Ij 


2j 



(16) 


Otherwise, if the strata contains a Croup III element that is 
dependent, upon the Croup I and Croup II components of the strata, 
the Croup III component is not individually computed, but shall 
be accounted for in the total strata variance computations in 
paragraph 3.7. .2.4. 

3. 2. 2.4 Ratio Area Est imate V ari ance . There are three categories 
into which the Group II i data may fail, each producing a uniquely 
determined variance estimate for the strata. The first category 
shall include that Croup III data which was estimated from the 
Group I and Croup II elements within the strata. The second 
category shall include that. Croup III data which was estimated 
fron all the Croup I and Group !I data within some space (zone, 
region, or country) as well as the Croup I and Group ii data 
within the strata. The third category shall consist of Group III 
data that was estimated from all the Group I and Group II data 
within some space (s) but no Group I or Group II data within the 
strata . 


A. Category 1. Within Strata Variance: 


and 


V(A. J ) 


V 

' j 


V 


(V, .♦!' 

1 ) 2.1 




2 <V A 


2 j 


(17) 


Where : 


V 

V 


lj 

2 j 


Group I area estimate variance for the jtk strata 
Group ! I area estimate variance for the j th strata 


3-15 



JSC-10009 


A^ = Group I area estimate for the jiTz strata 

A = Group II area estimate for the jth strata 
• 2} 

x. = (1/m) (r^+l) 

J V=1 

HI . 7 

. 2 =-[l/(m-l)3 £ (r 3 -r J )^ if m > 2 

CT. V=1 

1 

= 0 if m = 1 

m = number of years available 

r ^ = [ZW ^/(SW-^+SW ^) ] for the vth year 
j X z 

m 

r 3 = (1/m) V r 3 
v=l 

The superscripts imply "computed from data belonging to 
that space." 

B. Category 2. Between Space Variance : 

V(A'. S ) = V. 

3 .3 

and 

V. = (V^+V^) + 2(A x S +A 2 S ) 2 (<t S ) 2 

+ (V 1 S +V 2 S )(t S ) 2 + 2(V lj+ V 2j ) r s (18) 
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Where : 


V X S =■ V(Aj S ) 

V - vca 2 s ) 


(t S ) 2 


(1/n) 




2 


o . 
J 


2 


3 [1 / (m- 1 ) ] ^ ( r^ S - r s ) 2 if n > 2 
- 0 i f m ■ 1 


f s = (1/m) V r v S 

■3-1 


The space (s) includes all data within the level from which 
Group III estimate was obtained (i.c., ’.one, region, country). 

C . Ca t eg o ry 3 . Withon t S t r a t a Vari a n_c_e : 


and 


V(A. S ) 



V. = (VjMyu 


S)2 


2(A 1 S *A : S )(c S ) 2 


(19) 


The space (s) includes all data within the level from 
which the Group III estimate was obtained (i.c., zone, 
region, country) excluding the strata in question. 


3-15 



JSC-10009 


3.2.3 Aggrcg.i t ion C ontrol Logic . Spring and winter wheat area 
estimates at the strata level shall be aggregated separately to 
the zone, region, or country level based upon user request. The 
associated variance of the area estimates shall be aggregated 
along with’ each area estimate. 


The total wheat area estimate and the associated variance shall 
be computed as the sum of the spring and winter wheat components 
of the estimates for a given crop year (CY). 


a totai. , 

CY " A SPRING, 

CY 

A W INTER, CY 

(20) 

V TOTAt. , 

CY " V SPRING, 

CY 

'winter, CY 

(21) 

The standard statist 
shall be computed as 

ics associated with a total wheat 
defined in paragraph 3.2.2 from 

area estimate 
the total 


variance estimate. 


3 . 2 . 3 . 1 Area A ggregat ion to the Zo ne, Region, and Co untry L evels 
for a Sprin g/ W inte r ( : r op . Spring a nd w i n r. e r w h cat area e ? t i m a t e s 
for a tone level element shall be obtained through aggregation 
of the area strata within the tone element of the appropriate crop 
type. The aggregation follows similarly for each element in the 
region and country levels. The computation used for the area 
estimate aggregation shall be as follows for both spring and win- 
ter wheat crop. 


A . Area lis t in at e at Z on e L eve l : 


Where : 


n 


2 


A u 
nmf. 



A . . 

nm j- i 


( 22 ) 


Number of strata in the Hh zone, mth region, 
nth country 


A a • 

nmHj 


Area estimate of the jth strata, tin zone, nth 
region, nth country. 
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B. 


Area Estimate a t Reg i on level : 


A 

nm 



^nnt 


(23 ) 


Where : 


n R = Number of zones in the m*J: region, nt./: country. 
C . Area bs_r j natc at Coun t rv Le ve 1 : 


A 

n 



m- 1 


(24) 


Where : 


r,^ - Number of regions in the n r h country. 

The area aggregation for the spring and winter wheat crops shall 
be performed separately. 

3.2. 3. 2 Va r i a nc e A g g r c t; a t i on . The estimate of error (variance) 
associated with each area estimate at the strata level shall be 
aggregated in parallel with each area estimate aggregation so that 
there is a variance associated with each estimate at all hierarchi 
cal levels. The associated variance shall be used to compute the 
standard statistics as defined in paragraph 3.2.2 at each hier- 
archical level. 


A. 


V ar i ance A g g r e g a t i or._to th e Zone Lc v e 1 . The v a r i a n c e 
of the strata area estimate shall be aggregated to the 
zone level in one of two ways depending upon the existence 
and source of Group III data. If there is at least one 
strata within the zone in question that has a Croup III 
element which was estimated from region or country data, 
then the variance of the area estimate of the zone shall 
he calculated as follows: 


V (A 


nnu 


- V 

nmZ 
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and 


V c « V . 
. nmx nmx- 


si 


M 




Nhere : 
si 


s 2 


s 3 


V 


si 


nmi 


4 2 (A . „ s2 +A, 

1 ,t 2 , t 


S V (o'/ 5 ) 2 ] 


(25) 


M 

V 

V 


( t / 3 ) 2 


Space containing all the Group I anil Group II 
elements wit a in the it/: tone and all Group III 
elements that were not estimated hv a region 
or country ratio (i.e., 'reduced -one') 

r Space containing all strata within the tth zon 
except the tth strata (i.e., that strata being 
estimated) 

= Space containing all Group I and Group I! ele- 
ments not in t he hr h zone hut within the level 
(nrh region or nth country; iron which the 
Group III element ir. the tth strata was esti- 
ma t ed 

" Variance of the 'reduced zone’ (si) of the 
tth zone, mth region, nth country 

J Number of Group III elements within trie /.th- 
rone that were estimated from data outside 


s2 

tb.e tth 

zone 

l.t 

n Variance 

c f 

s 2 

2 . t 

- Variance 

o f t 


l.t 


s 2 


2 1 1 


rt 


W) E (r v .‘V 

v^l » v 


m « Number of years of historical data avail- 
able 
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p [K ^ J (W ^ S ^ ) ) ^ or vt ' 1 year 

Group ! area estimates frun the space s2 (docs 
not contain the tth strata) 


t = Group II area estimates from the space s2 (does 
* not contain the tih strata) 




F ■ System constant scale factor (default " I) 

HI - Space containing Group 111 elements estimated 
from within t he ith 7 . one 

1:2 - Space containing Group III elements estimated 
from without the i :.h z one 
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Then 


V 


nm<t 


si 



V 


si 


Rinfc j 


(26) 


2 . 


Where : 


V 


si 


nraij 


Variance of the area estimate in the jth 
strata of the space si in the it/: zone, 
nth region, nth country 


N1 


If 


Number of strata in the ’reduced zone” 
(space si). 


E *:, 1 < F £ w,< 

t e E 2 t r. H 2 


Then 


V 


nmi 


si si si w si . 2 

k l.nnJ. 2 ,nm£ M ; 


Where : 

v 

V 1 ,nm£ 


si 


2 , nmZ 


2 (A 


si 


1 ,nm£ 


■A, 


,nmi 


S V 


s i 2 
(o S ') Z 


(27) 


Group 1 variance of the 'reduced zone' 
(space si) tik zone, nth region, nth 

Country- 

Similar Group II variance 
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si 


l,nm& 


si 


l 2 ,nm£ 


Group I area estimate of the ’reduced 
zone' (space si) of the Ith zone, m th 
region, nth country 

Similar Group II area estimate 


(t 



2 


, si. 2 

O ) 


m . ? 

- U/nO Y. (r, *1) 

V=1 

m it? 

= [1/ (m-1) ] 52 ( r v S ' ?S ) if m > 2 

v— 1 

= 0 if m = 1 

sl si si si 

r v " 0*3 )/(W l +W 2 J 

f sl - Cl/.) t \ Sl 

V=1 


Otherwise, if there are no Group III elements within the 
zone that were estimated based upon data from outside the 
zone, the variance of the area estimate of the zone shall 
be calculated as follows : 


N2 


nmJl 


- E 

:=i 


V 


nm$,j 


(28) 


Where : 

V . = Variance of the area estimate of the j th strata, 
^ %th zone, m th region, nth country 

N2 = Number of strata in the Ith zone, m th region, 
nth country. 
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B . Variance Aggrcgat i or^ to the Reg ion and Cc u r. t ry Level s . 

The variance aggregation to the region and country levels 
shall he performed as follows: 


1 . 


For variance 


at the region level: 


'R 

v - y* v „ 

nm nmfc 


(29) 


2. For variance at the country level: 


V 

n 


r c 

■ E v 

m- 1 


(30) 


The between zone variance component shall be considered 
as negligible. 

3. 2. 3. 3 Proc ess Spe cif ication Par amete rs. Program parameters 
shall be e n"t e red wh" i c h Vi 1 1 govern the processing performed and 
the outputs generated. The following subsections define the pro- 
gram parameters which r.ay be entered interactively or in batch 
node . 

A. Cou nt r y . One of the LACIE countries must be selected. 

B. Spr i ng/ Wintcr/Tota 1 . A designation of the crop type 
(spr i ng/wi nt. e r/ 1 o t a 1 ) shall he required tc identify the 
desired processing. 

C. Cr op Year. The crop year of data to be processed must be 
spec i f ied. 

i). Level/ HI cnc r.t . Both the hierarchical level (strata, zone, 
region, country) to which aggregation is to be performed 
and the specific element within that level for which 
aggregation is to be perferned mist be specified. 

H. Classif ica t. ion Res ult s Criteria. A threshold level of 
classification results may he entered (range 0-15). 
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F. Ind i v i d ua 1 /Group R eport . Output of results nay be an 
individual report for the highest level of aggregation or 
a group of reports for the highest level of aggregation 
and all interim aggregations. 

G. Data Sup pression Flag . Optional suppression of the 
supportive and historical report data fron being output 
may be requested. If no specification is made, the de- 
fault. of no suppression shall he assumed. 

H. Output D evice's). The CAS line printer, CRT, or a magnetic 
tape may be designated as an output device for report pur- 
poses . 

I. Biostag c T hr esh old. This entry shall define the upper 
level and lower level of the biological stage (s) of the 
crop within a crop year for the data to be processed. 

J. Acqu i s i t. i on Date Th resh old. This entry shall define the 
upper level and lower level of the acquisition dale(s) 
within a crop year of the LAND SAT data to be processed. 

3.2.4 Input Requi remen ts. Input data necessary to support the 
area estimation, aggregation, and reporting requirements shall 
consist of historical data, classification data, and allocation 
data. 


3.2.4 


1 C lassif i cat ion D ata . The following data shall be requ 
Estimated wheat area per sample segment in hectares 


Sample segment identification number 
Percentage classified as wheat 
Percentage classified as other 


Percentage classified as unidentifiable 
Percentage threshol ded 


red 
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• Identification as spring or winter wheat sample segment 

• Biological stage of each acquisition 

• Date of LANDSAT acquisition (s) 

• Date of classification 

• Classification evaluation code 

• Crop year 

• Number of acquisitions used in classification. 

3. 2.4. 2 Historical Data . The following historical data shall be 

required: Spring and winter historical wheat area for the last 

15 years for each substrata/strata, zone, region, and country. 

One of the IS years must be the epoch year for each country. 

3. 2.4. 3 Allocation Data. The CAS data base shall require the 
following data to support the CAS estimation and aggregation 
formulations . 

A. For the lowest, hierarchical element of each country: 

• Identification number 

• Number of agricultural segments 

• Nominal allocation of element • Group i, 11, or III 

• Number and identification of sample segments allocated 
to the element 

• Ratio of true element, area to gross pseudo- e 1 crr.cn t area 

• Indicator cf whether spring, winter, or both spring and 
winter wheat is grown in the clement. 
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B. For each I.AC Ii: country: 

• Identification ar.d correlation of hierarchical elements 
. (regions, cones, strata, substrata) 

• Indicator of density of historical data (bilevel) 

• Identification of epoch year data (data used in the 
allocation process). 

3 . 2 . S Outpu t Rcqu i_r_(?nu: n t. s . Th e vh eat a r e a 
statistics at the base level and succeeding 
area aggregation shall be stored for future 
to the output device (s). 

3. 2. 5.1 Report Generation. The area estimates and standard 
statistics on the strata, cone, region, and country levels shall 
be available for presentation in conjunction with historical data 
for all elements within the requested aggregation level. These 
reports may be generated or. a batch device and/or CRT terminal. 
The contents for area estimate reports are specified in paragraph 
3.5. In addition, auxiliary reports may he generated containing 
substrata area estimate information. 


estimates and standard 
levels as a result of 
reference and reported 


3 . 1 , 5 . 2 Res ul t s _S t o r a g e . Interim s t o r a g e <• f w h e a 
and associated variances at the strata level shall 
use by the wheat production processing phase. 


area estimates 
be done for 
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3.3 YIELD ESTIMATION 

Yield estimates, along with their associated variances, shall be 
given at the strata level for all strata. The yield estimate for 
any higher element in the hierarchy shall be computed from the 
area and production estimates for that element. These yield esti- 
mates, along with the associated historical and statistical infor- 
mation, shall be output as requested in standard report formats. 

The input data- shall consist of wheat area and production esti- 
mates, yield data, historical data, and allocation data. The 
output reports shall contain yield estimates with statistics as 
specified in paragraph 3.3.5. 

3.3.1' Yield Estimates . Wheat yield estimation shall be computed 
at the zone, region, and country levels for each LACIE country 
as requested. Separate estimates shall be performed for the 
spring and winter wheat crop, while no total wheat yield estimate 
is produced. 

3. 3. 1.1 Yield Estimates at Strata Level- . Yield estimates shall 
not be computed at the strata level since these data are provided 
monthly by the Yield Estimation Subsystem (YES). These estimates 
shall be used in production estimate and variance computations on 
the strata level (paragraphs 3. 4. 1.1 and 3.4. 2.1) and shall be 
reported unchanged. 

3.3. 1.2 Yield Estimates at Zone, Region, and Country Levels . The 
yield estimates for each element (e) within the zone, region, and 
country levels shall be computed as follows : 

Y = P /A (31) 

e e e 


Where : 

e = Any zone, region, or country 

A - Area estimate for element e 
e 

P = Production estimate for element e (as defined in paragraph 
6 3.4). 
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3.3.2 S tand ard S tatist ics. The standard statistics tor yield 
estimation are composed of standard deviation (0. coefficient of 
variation (CV) , 90 percent confidence limits (Cl.) > and probability 
of achieving less than 10 percent error (»'l) . Each statistic is 
a function cf variance V as follows: 

o s = JT (32) 

Where s - Any area clement. 

CV - a / Y (33) 

s s s 

Where Y = Yield estimate of the sth element, 
s 

CL (+ ' - Y ♦ (1.645) o 
s s s 

CL s < -’- 1 - Y s - (1.64S) a (34) 

= erf j(l/ N ")(l/(lCCV s )lJ (35) 

Where : 

erf (8/\ (2) - Error function associated with the normal prob- 
ability curve 

0 ■ Standardised random variable. 

3. 3. 2.1 Yield E stim ate Variance at Strata Le vel . The variance 
of the yield estimates at the strata level are not computed since 
these data shall be provided monthly by the Yield Estimation Sub- 
system (YES). These estimates shall be used in production vari- 
ance computations on the strata level (paragraph 3. 4. 2.1) and are 
reported unchanged. 
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3. 3. 2. 2 Yield Estimate Variance at Zone, Region, and Country 
Levels . The variance estimate associated with each element (e) 
above the strata level is computed as follows: 


voy 


e 

L 


E V <V E V <V E [ Y ej v <v] 



e 


(36) 


•Where : 
Y 

e 

• L 

P . 
ej 


V(P .) 


ej 

V (A . ) 
ej 


ej 


= Yield estimate for element e 

= Number of strata (j ) within the area element e 

= Production estimate of the j th strata in element e 

- Variance of P . (as defined in paragraph 3.4. 2.1) 

ej 

= Area estimate of the j th strata in element e 

= Variance of A . 

ej 

= Yield estimate of the j th strata in element e 


3.3.3 Process Specification Parameters . In order to govern the 
processing performed and the outputs generated, certain program 
parameters shall be entered. The following paragraphs define 
those program parameters. 


A. Country . One of the LACIE countries must be selected. 

B. Spring/Winter . A designation of the crop type (spring/ 
winter) shall be required to identify the desired proc- 
essing . 

C* Crop Year . The crop year of data to be processed must be 
specified. 
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D. Lcvc 1 / Elem ent . Both the hierarchical level to which proc- 
essing is to be performed and the specific element to 
process within that level must be specified. 

K. Classification Re sts I ts Criteri a.* A threshold level of 
classification results may be entered (range 0-15). 

F. Ind iv idua 1 /Crou p Report . Output of results may be an 
individual report for the highest level processed or a 
group of reports for the highest level processed and all 
subclenents of the requested element. 

G. Data Suppression Fl ag . Optional suppression of the sup- 
portive and historical report data from being output may 
be requested. If no specification is made, the default 
of no suppression shall be assumed. 

H. Output De vice (s). A CRT and/or hatch device may he des- 
ignated as an output device for report purposes. 

I. Biost agc Th res ho I d.* This entry shall define the upper 
level and lower level of the biological stage (s) of the 
crop within a crop year for the data to he processed. 

J. A cquisiti on Date Th resho ld.* This entry shall define the 
upper level and lower level of the acquisition datc(s) 
within a crop year of the LAND SAT data to be processed. 

K . Yie 1 d_ Fs t imat e_ Da t e Thres h oil. T h is entry shall d e f i r. e 
the upper and lower level date, wit); in a crop year, of 
the yield estimate data at the strata level to he used 
in processing. 

3.3.4 Input Rcu u i re cent s . Input data necessary to support the 
yield estimation and reporting requirements shall consist of 
yield data, historical data, allocation data, ar.d area ar.d pro- 
duction estimate data. 


*Thesc parameters shall be checked 
area and production estimate data 
requested for the yield estimates 


to assure that the associated 
arc of the same type as that. 
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3. 3. 4,1 Yield Data . The following data shall be required at 
each area strata: 

• Strata identification number 

• Spring/winter wheat crop indicator 

• Acquisition date 

• Crop year 

• Yield estimate value for specified wheat crop 

• Variance of yield estimate for specified wheat crop. 

.3. 3. 4. 2 Historical Data . Spring and winter historical wheat 
yield for the last 15 years for each strata, zone, region, and 
country shall be required. 

i 

3.3.4. 3 Allocation Data . Identification and correlation of 
hierarchical elements (regions, zones, strata) for each LACIE 
country shall be required. 

3. 3. 4. 4 Area Production Estimate Data . The following data shall 
be required for every zone, region, and country -element: 

• Area estimate 

• Variance of area estimate 

• Production estimate 

• Variance of production estimate. 

These area and production estimates shall have been specified in 
the same manner as that requested for the yield estimate (para- 
graph 3.3.3). 

3.3.5 Output Requirements . The wheat yield estimates and stand- 
ard statistics shall be stored for future reference and reported 
per request to the output device (s). These estimates shall be 
presented at the strata level and succeeding levels as a result 
of yield data processing. 
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3. 3. 5. 1 Report Gene ration . The yield estimates and standard 
statistics shall he available for presentation in conjunction with 
historical data for ail elements below the requested processing 
level. These reports may be generated on a magnetic tape, line 
printer, and/or CRT terminal. The contents of the yield estimate 
reports are specified in paragraph 5.5. 

5. 5. 5. 2 Results Sto rage . Interim storage of wheat yield esti- 
mates and associated variances at the strata level shall be done 
for use by the wheat production processing phase. 
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3.4 PRODUCTION ESTIMATION AND AGGREGATION 

The production applications software shall calculate production 
estimates and standard statistics from the input data (paragraph 
3.4.4). The production estimates and variances shall be aggre- 
gated to the level specified by the process specification param- 
eters (paragraph 3. 4. 3. 3). The production and statistical informa- 
tion shall be output (paragraph 3.4.5) as requested in standard 
report formats. The input data shall consist of area and yield 
estimate data, historical data, and allocation data. 

3.4.1 Production Estimates . Wheat production estimation shall be 
calculated at the strata level from area and yield estimates. 
Separate estimates shall be performed for spring and winter wheat. 

Production estimation shall be performed on a base level of strata 
elements. The production estimate shall be calculated by the 
following, equation. 


Where : 


P. = A. Y. 
J 3 1 


(37) 


A. = Wheat area in the j th strata 
3 

Yj = Wheat yield in the j th strata 

3.4.2 Standard Statistics . The standard statistics for pro- 
duction estimation and aggregation are composed of standard 
deviation (o) , coefficient of variation (CV) , 90 percent con- 
fidence limits (CL), and probability of <_ 10 percent error (ft) . 
Each statistic is a function of variance (V) as follows: 


a 

s 



Where s = Any hierarchical element 


Where P 

s 


CV 


s 


a /P 
s 


s 


Production estimate of the sth element 


(38) 


(39) 
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- P g ♦ (1.645) o s 


= P s - (.1.645) o g 


i s - erf j(l/vT)[l/(10CV s )J 


(41) 


erf (0/ v T) 


Hrror function associated with the normal prob- 
ability curve 


0 = Standardized random variable 

3 . 4 . 2 . 1 Str ata I. eve 1 Pr oduc lion Kst in_a te _Var_i_ar.ee. The var i ance 

of wheat production estimates at the strata level is giver, by: 

X. - V(P.) - V(A ) (V(Y.) ♦ Y. 2 ] ♦ [V(Y ) A. 2 ] (42) 

Where : 


V(Aj) " Variance of wheat area in the j th strata 
V(Yj) ■ Variance of wheat yield in the j in strata 

A. and Y. are as previously defined. 

J ) * 

3.4.3 Ark r eg a t i on Con trol i.ogic. Spring and winter wheat pro 
duction estimates at the strata level shall he aggregated sep- 
arately to the zone, region, or country level based upon user 
request. The associated variance of the produc; ion estimates 
shall he aggregated along with each production estimate. 


The total wheat production estimate and the associated variance 
shall he computed as the sum of the spring and winter wheat com- 
ponents of the estimates for a given crop year (CY). 
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P TOTAI. , CY 

^SPRING , CY * P W1NTF.R, 

CY 

(43) 

V TOTAL, CY 

^ SPRING , CY + V W ! INTER , 

CY 

(44) 


The standard statistics associated with a total wheat production 
estimate shall be computed as defined in paragraph 3.4,2 from 
the total variance estimate. 

3 .4.3.1 Product ion Agg rc gaf i o n t o the Zo ne , R eg ion , and (tonr. tr/ 
Leve ls ♦ Spring and winter wheat production estimates for a zone 
level element shall be obtained through aggregation of the pro- 
duction strata within the zone element of the appropriate crop 
type. The aggregation shall follow similarly for each element 
in the region and country levels. The computation used for the 
production estimate aggregation shall be as follows for both spring 
and winter wheat crops. 

A . Produ c t_i o n list i rate at the Zone L evel : 


N 


nnu. 1 


- E 


nm ?. j 


(45) 


j-1 


Wh ere: 


N 


Number of strata in the lih zone, mth region, nr/; 
count ry 




nmtj 


Production estimate of the j th strata, Ith zone, 
mt/i region, nth country 


B. Production Hstimate at the Region Level: 


n r 

r . V p . 

nm ' nra 

1 


(46) 


Where : 


N - Number of zones in the mfh region, nth country 
K 
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C. 


Production Kstimat c at the Country Level : 


P 

n 



P 

nn 


(47) 


Where : 


N c ° Number of regions in the nth country 

The production aggregation for the spring and winter wheat crops 
shall be perforned separately. 

3.4. 3.2 Var ia nce Aggregat io n to the Zone, He gion , and Co untry 
Lcvc 1 s . The estimate of error (variance) associated with each 
production estimate at the strata level shall he aggregated in 
parallel with each production estimate aggregation so that there 
is a variance associated with each estimate at all hierarchical 
levels. The associated variance shall be used to conputc the 
standard statistics as defined in paragraph 3.4.2 at each hier- 
archical level. 


The variance aggregation tc the zone , region, and country levels 
shall be performed as fellows for both spring and winter wheat 
crops . 

A. Production Variance at the Zone Level: 


S„ 

V - V(P 1 - V V „ . (4 3) 

nn£ nm t nnt; v ' 

j-1 

Where : 

N = Number of strata in the ilk zone, mi*: region, 

nth country 

V fl . * V(P__.i) Variance of the production estimate of 
nn i j . '"I'R ) . , . , 

the jr /: strata, t-tn zone, mi/< region, nth 

country 
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Where : 

N - Number of regions in the nt/: country 

3 . 4 . 3 . 3 Pr ocos s Spec i f i ca t ion Par anet er_s . The foil owing para - 
graphs define the process specification parameters which may ho 
entered interactively or in hatch node. The process specification 
parameters shall he used to define the processing that is required 
and the output products that are generated. 

A. C ount ry . One of the BACH: countries must be selected. 

B. S pr i ng / Wi n t e r /T o t a 1 . A designation of the crop type 
(spr ing/winier/total) shall be required to identify the 
desired processing. 

C. C rop Year. The crop year of data to be processed must be 
specified . 

I). l.eve l/lilencnt . Both the hierarchical level (strata, ?.cr.c , 
region, country'; tc which aggregation is to be performed 
and the specific clement within that level lor which 
aggregation is to he performed rust he specified. 
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E. Classification Results Criter ia.* A threshold level of 
classification results nay be entered (range 0-15). 

F. Individua l /G roup Report . Output of results may be an 
individual report ter the highest level of aggregation or 
a group of reports for the highest level of aggregation 
and all interim aggregations. 

G. Data Suppress i on Flag . Optional suppression of the sup- 
portive and historical report data from being output may 
be requested. If no specification is made, the default 
of r.o suppression shall be assumed. 

H. Output Device (s). A CRT and/or a batch device nay be 
designated as an output device for report purposes. 

I. Bios tage Thr csho 1 d . * This entry shall define the upper 
level and lower level of the biological stage(s) of the 
crop within a crop year for the data to be processed. 

J . Ac quis i t i_on i) at c_ Thres hold . * This entry shall define t he 
upper level and lower level of the acquisition date(s) 
within a crop year of the I.ANDSAT data to bo processed. 

K. Yield Es timate Date Threshold.* This entry shall define 
the upper and lower level date, within a crop year, of 
yield estimate data at the strata level to be used in 
processing. 

3.4.4 Input Re quirem ents. Input data necessary to support the 
production estimation and reporting requirements shall consist 
of wheat area data, yield data, historical data, and allocation 
data. 

3.4.4. 1 Area and Yi eld Estima t e D ata. Wheat area and yield data 
shall be required to compute the production estimates and variances. 

*These parameters shall be checked to assure that the associated 
area and yield estimate data are of the same type as that re- 
quested for the production estimates. 
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The area and yield data shall be checked for compliance with the 
process specification parameters listed in paragraph 3. 4. 3. 3. The 
follovring data shall be required: 

• Estimated area per strata 

• Estimated yield per strata 

• Area estimate variance per strata 

• Yield estimate variance per strata. 

3. 4. 4. 2 Historical Data . Spring and winter historical wheat 
production data for the last 15 years for each strata, zone, 
•region, and country shall be required for report purposes. 

3; 4. 4. 3 Allocation Data . Production estimation and aggregation 
shall require data which identifies each LACIE country's hierar- 
chical elements (region, zone, strata). 

3.4.5 Output Requirements . The wheat production estimates and 
standard statistics at the base level and succeeding levels as 
a result of production aggregation shall be stored for future 
reference and reported to the output device (s). 

3.4. 5.1 Report Generation . The production estimates and stand- 
ard statistics on the strata, zone, region, and country levels 
shall be available for presentation in conjunction with histor- 
ical data for all elements within the requested aggregation level. 
These reports may be generated on a batch device and/or CRT ter- 
minal. The contents of -the production estimate reports are 
specified in paragraph 3.5 

3.4. 5. 2 Results Storage . Interim storage of wheat production 
estimates and associated variances at each level shall be done 
for use by the wheat yield processing phase. 
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3 . S REPORT DATA RF.QU I REM ENTS 

Reports which were generated through the area, yield, or pro- 
duction estimate processing phase may be displayed on the request 
cd output devicc(s). The »avai lablc devices are CRT, magnetic 
tape, and line printer as defined in the System Software Section 
(section 2). The estimation reports arc presented in cither 
metric or USA units based upon user request at input command time 
(paragraph 2.2.1). Following is a cross-reference of the unit 
specifications that are presented for each data type. Separate 
reports shall be generated for spring, winter, and total crop 
types . 


Data Ty pe 


Metric 


USA 


Area 

Yield 

Production 


Hectares Acres 

Quintal s/Hcctarc Bushels/Acrc 

Thousand Metric Tons Thousand Bushels 


Additional reports containing yield strata data or Worldwide 
Meteoriogical Organization (WMO) position data may be output per 
a request at input command time (paragraph 2.2.1). All report 
data that is available in the data base shall be presented in 
each of the report requests. 

3.5.1 Spring/ Win ter /Tot a_l_ >.'h eat Area hs t i mat i on Rep or t_s . The 
area estimation reports shall be output as strata, zone, region, 
or country reports. The parameters required for each report are 
specified in the following paragraphs. 

3. 5. 1.1 S trata . Each strata area estimation report shall consis 
of the following information: 

A. Head ing 

1. Strata number 

2. Crop year 

3. Units of measure 
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4. Thresholds applied on the report data 

• Minimum classification evaluation code 

• Biological stage range 

• Acquisition date range 

5. Associated zone, region, country 

6. Title 

7. Spring cr winter or total. 

B . S ecurity B! o c_k 

1. CAS monthly report date - user input 

2. Date through which maximum protection must be retained 

3. Date through which restricted access must he retained 

4. Date of report generation. 

C . Are a Estimate for Strata 

1. Area estimate of strata 

2. Estimated standard deviation 

3. Probability of achieving 10 percent or less error 

4. Coefficient of variation 

5. 90 percent confidence limits. 

D . S uppor t ivc j ) a t. a * 

1. Number of sample segments allocated 

2. Number of sample segments used in estimate 

3. ID's of sample segments used in estimate. 

^Output if base hierarchical level is strata. 
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H . Supportive Dat a* 

1. Number of Croup I substrata and ID's 

2. Number of Grou*p II substrata and ID's • 

3. Group I estimate and standard deviation 

4. Group II estimate and standard deviation 

5. Group III estimate. 

F. Historical Ar ea D ata 

1. Number of years of data 

2. Maximum area and year 

3. Minimum area and year 

4. Mean area 

5. Standard deviation. 

G. Comme ntar y - 160 characters of commentary narrative. 

3.S.1.2 Zone. The zone area estimation report shall consist 
of the following information: 

A. He ading 

1. Zone number 

2. Crop year 

3. Units of measure 

4. Thresholds applied on the report data 

• Minimum classification evaluation code 

• Biological stage range 

• Acquisition date range 

^Output if base hierarchical level is substrata. 
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5. Associated region and country 

6 . Title 

7/ Spring or winter or total. 

B. Security Bloc k 

1. CAS monthly report date - user input 

2. Date through which maximum protection must he retained 

3. Date through which restricted access must be retained 

4. Date of report generation. 

C. Area Hst imafe f or Zone 

1. Area estimate of zone 

2. Hstimated standard deviation 

3. Probability of achieving 1C percent or less error 

4. Coefficient of variation 

5. 90 percent confidence limits. 

D. Support i ve Data (for each stratum in the zone) 

1. Strata ID 

2. Area of strata 

3. Standard deviation of strata 
K. Histo rical Data 

1. Number of years of data 

2. Maximum area and year 

3. Minimum area and year 
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4. Mean area 

5. Standard deviation. 

F. Comr.cn t ary - 160 .characters of commentary narrative 

3. 5. 1.3 Region . Each region area estimation report shall 
sist of the following information: 

A. Heading 

1. Region name 

2. Crop year 

3. Units of measure 

4. Thresholds applied on the report data 

• Minimum classification evaluation code 

• Biological stage range 

• Acquisition date range 

5. Country 

6. Title 

7. Spring or winter or total. 

B . Securi ty B lock 

1. CAS monthly report date - user input 

2. Date through which maximum protection must he r 

3. Date through which restricted access must be re 

4. Date of report generation. 


con- 


eta ined 
t a i ned 
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C . Area Kstisnate for Regio n 

1. Area estimate of region 

2. Hstimated standard deviation 

3. Probability of a 10 percent or less error 

4. Coefficient of variation 

5. 90 percent confidence limits. 

D. S upportive Data (for each zone in the region) 

1. Zone ID 

2. Area of tone 

3. Standard deviation of tone 

0 . Histo rical Dat a 

1. Number of years of data 

2. Maximum area and year 

3. Minimum area and year 

4. Mean area 

5. Standard deviation 

F. Conne nt_arv - 160 characters of commentary narrative. 

3.3.1. 4 Country. each country area estimation report shall con- 
sist of the foil cw : ng i r. f orma t i on : 

A . Head i ng 

1. Country name 

2. Crop year 
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3. Units of measure 

4. Thresholds applied or, the report data 

• Minimum classification evaluation code 

• Biological stage range 

• Acquisition date range 

5. Title 

6. Spring or winter or total. 

B. Securit y Bloc_h 

1. CAS monthly report date - user input 

2. Date through which maximum protection must be retained 

3. Date through which restricted access must be retained 

4. Date of report generation. 

C. Area list i mate for C oun try 

1. Area estimate of country 

2. Estimated standard deviation 

3. Probability of achieving 10 percent or less error 

4 . Coefficient of variation 

5. 90 percent confidcr.ce limits. 

I). Supportive Data (for each region in the country) 

1. Region ID 

2. Area of region 

3. Standard deviation of region. 
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F.. Historical Da ta 

1. Number of years of data 

2. Maximum area and year 

3. Minimum area and year 

4. Mean area 

5. Standard deviation. 

F. Comme ntary - 160 characters of commentary narrative. 

3 . S . 2 Spr i ri£ / W i n t or Wheat V i eld F s lim a t i on Repo rt s . The yield 
estimation results shall be output as separate spring wheat or 
winter wheat crop reports at the strata, zone, region, or country 
levels. The parameters to be presented in each report are spec- 
ified in the following paragraphs. 

3. 5. 2.1 Str ata. each strata yield estimate report shall consist 
of the following information: 

A . He a d ing 

1. Title 

2. Strata number 

3. Strata name 

4. Zone, region, and country to which the strata belongs 

6. Crop year 

6 . Crop type (spring or winter) 

7 . Units indicator 

8. Thresholds applied on the report data 

• Minimum classification evaluation code 

• Biological stage range 
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• Acquisition date range 

• Yield estimate acquisition date range 

• V 

R. Security Block 

* - — 

1. CAS monthly report date - user input 

2. Date through which maximum protection must be retained 

3. Date through which restricted access must be retained 

4. Date of report generation. 

C . Yield Estina t e f o r Str a t a 

1. Strata yield estimate 

2. Standard deviation of yield estimate 

3. Coefficient of variation of yield estimate 

4. Upper and lower confidence limits c: yield estimate 

5. Probability of achieving less than 10 percent, error 
in yield estimate. 

D. Historical Yi eld Data 

1. Number of years of historical data available 

2. Maximum historical yield for the strata with the associ- 
ated year 

3. Minimum historical yield for the strata with the associ- 
ated year 

4. Mean historical yield for the strata 

5. Standard deviation of historical yield. 

H. Commentary - 160 characters of user input commentary. 
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3. 5. 2. 2 Zone. F.ach zone yield estimate report shall consist of 
the following information: 

A . Head ing 

1. Title 

2. Zone number 

3. Zone nanc 

4. . Region and country to which the zone belongs 

5. Crop year 

6. Crop type (spring or winter) 

7. Units indicator 

8. Thresholds applied on the report data 

• Minimum classification evaluation code 

• Biological stage range 

• Acquisition date range 

• Yield estimate acquisition date range 

B . S e cur it. y F> lock 

1. CAS monthly report date - user input 

2. Date through which maximum protection must be retained 

3. Date through which restricted access must be retained 

4. Date of report generation. 

C . Yie ld E s t i ma te f o r Zone 

1. Zone yield estimate 

2. Standard deviation of yield estimate 
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3. Coefficient of variation of yield estimate 

4. Upper and lower confidence limits of yield estimate 

5. Probability of achieving less than 10 percent error 
in yield e s t ima t e . 

I) . Support : v c_D a t a 

1. Zone area estimate 

2. Zone area standard deviation 

3. Date of zone area estimate calculation 

4. Zone production estimate 

5. Zone production standard deviation 

6. Date of zone production estimate calculation. 

P . Hi s tor i c al Yield D ata 

1. Number of years of historical data available 

2. Maximum historical yield for the zone with the associ 
a ted year 

3. Minimum historical yield for the zone with the associ 
a ted year 

4. Mean historical yield for the zone 

5. Standard deviation cf historical yield. 

F. C omme ntary - 160 characters of user input commentary. 

3. 5.2. 3 Rcgi or. . bach region yield estimate report shall consis 
of the following information: 

A . H ead j ng 

1. Title 

2. Region number 
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3. Region nano 

4. Country to which the region belongs 

5. Crop year 

6. Crop type (spring or winter) 

7. Units indicator 

8. Thresholds applied on the report data 

• Minimum classification evaluation code 

• Biological stage range 
t Acquisition date range 

• Yield estimate acquisition date range. 

B. Secur i ty Block 

1. CAS monthly report, date • user input 

2. Date through which maximum protection must be retained 

3. Date through which restricted access must be retained 

4. Date of report generation. 

C . Yield Hs timat e for Re gion 

1. Region yield estimate 

2. Standard deviation of yield estimate 

3. Coefficient of variation of yield estimate 

4. Upper and lower confidence limits of yield estimate 

5. Probability of achieving less than 10 percent error in 
yield estimate. 
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D. Supportive Data 

1. Region area estimate 
2 : Region area standard deviation 

3. Date of region area estimate calculation 

4. Region production estimate 

5. Region production standard deviation 

6. Date of region production estimate calculation, 
t: . His tor i cal Yield D ata 

1. Number of years of historical data available 

2. Maximum historical yield for the region with the assoc 
a ted year 

3. Minimum historical yield for the region with the assoc 
ated year 

4. Mean historical yield for the region 

5. Standard deviation of historical yield. 

i ; . Comme ntary - 160 characters of user input commentary. 

3.5. 2.4 Count ry . bach LAC IK country yield estimate report shall 
consist of the following information: 

A . Heading. 

1. Title 

2. Country number 

3. Country name 

4 . Crop year 
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5. Crop type (spring or winter) 

6. Units indicator 

7. Thresholds applied on the report data 

• Minimum classification evaluation code 

• Biological stage range 

• Acquisition date range 

• Yield estimate acquisition date range. 

B. Security Block 

. ~ —A 

1. CAS monthly report date - user input 

2. Date through which maximum protection must be retained 

3. Date through which restricted access must be retained 

4. Date of report generation. 

C . Yield Hs t i mate f o r Count r y 

1. Country yield estimate 

2. Standard deviation of yield estimate 

3. Coefficient of variation of yield estimate 

4. Upper and lower confidence limits of yield estimate 

5. Probability of achieving iess than 10 percent error in 
yield estimate. 

D . Support iv e_ _Da t a 

1. Country area estimate 

2. Country area standard deviation 
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3. Date of country area estimate calculation 

4. Country production estimate 

5. Country production standard deviation' 

6. Date of country production estimate calculation. 

K . Historica l Yield Data 

1. Number of years of historical data available 

2. Maximum historical yield for the country with the 
associated year 

3. Minimum historical yield for the country with the 
associated year 

4. Mean historical yieid for the country 
3. Standard deviation of historical yield. 

F. Co mmen t ary - 160 characters of user input commentary. 

3.S.3 Spr i n^/lV i r.t or /To t a 1 Wheat Production Ms t i ma t i on Reports. 
The production estimate results shall be output as separate 
spring, winter, or total wheat, crop reports at the strata, zone, 
region, or country levels. The parameters to be presented in 
each report are specified in the following paragraphs. 

3.5.3. 1 S trata . Mach strata production estimate report shall 
consist of the following information: 

A . H ead i ng 


1. 

Title 


2. 

Strata 

number 

3 . 

Strata 

name 
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4. Zone, region, and country to which the strata belongs 

5. Crop year 

6. Crop type (spring, winter, or total) 

7. Units indicator 

8. Thresholds applied on the report data 

• Minimum classification evaluation code 

• Biological stage range 

• Acquisition date range 

• Yield estimate acquisition date range. 

B . Securit y Bloc k 

1. CAS monthly report date - user input 

2. Date through which rr.axinum protection must be retained 

3. Date through which restricted access must, be retained 

4. Date of report generation. 

C . Est i mate for S t rata 

1. Strata production estimate 

2. Standard deviation of production estimate 

3. Coefficient of variation o:' production estimate 

4. Upper and lower confidence limits of production estimate 

5. Probability of achieving less than 10 percent error or. 
production estimate. 
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D. Supportive Data 

1. Strata area estimate 

2 . Strata area standard deviation 

3. Date of area estimate calculation 

4. Strata yield estimate 

5. Strata yield standard deviation. 

E. Histori ca l Product io n Dat a 

1. Number of years of historical data available 

2. Maximum historical production for the strata with the 
associated year 

3. Minimum historical production for the strata with the 
associated year 

4. Mean historical production for the strata 

5. Standard deviation of historical production. 

F. Commen t ary - 160 characters of user input commentary. 

3. 5. 5. 2 Zone . Each zone production estimate report shall consist 
of the following information: 

A . He ading 

1. Title 

2. Zone number 
5. Zone name 

4. Region and country to which the tone belongs 
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5. Crop year 

6. Crop type (spring, winter, or total) 

7. Units indicator 

8. Thresholds applied on the report data 

• Minimum classification evaluation code 

• Biological stage range 
« Acquisition date range 

• Yield estimate acquisition date range. 

B. Security Block 

1. CAS monthly report date - user input 

2. Date through which maximum protection must be retained 

3. Date through which restricted access must be retained 

4. Date of report generation. 

C. Estimates for Zone 

1. Zone production estimate 

2. Standard deviation of production estimate 

3. Coefficient of variation of production estimate 

4. Upper and lower confidence limits of production esti- 
mate 

5. Probability of achieving less than 10 percent error in 
production estimate. 

D. Supportive Production Data (for each strata within the zone) 

1. Strata number 

2. Strata name 
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3. Strata production estimate 

4. Standard deviation of strata production estimate. 

E . Historical Produc t 'io n D a t a 

1. Number of years of historical data available 

2. Maximum historical production for the zone with the 
associated year 

3. Minimum historical production for the zone with the 
associated year 

4. Mean historical production for the zone 

5. Standard deviation of historical production. 

F. C ommen t ary - 160 characters of user input commentary. 

3. 5. 3. 3 Region. Each region production estimate report shall 
consist of the following information: 

A . Headin g 

1. Title 

2. Region number 

3. Region name 

4. Country to which the region belongs 

5. Crop year 

6. Crop type (spring, winter, or total) 

7. Units indicator 

S. Thresholds applied or. the report data 

• Minimum classification evaluation code 

• . Biological stage range 
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• Acquisition date range 

• Yield estimate acquisition date range. 

B . Security Bloc k 

1. CAS monthly report date - user input 

2. Date through which maximum protection must be retained 

3. Date through which restricted access must be retained 

4. Date of report generation. 

C . Est imates for Regi on 

1. Region production estimate 

2. Standard deviation of production estimate 

3. Coefficient of variation of production estimate 

4. Upper and lower confidence limits of production estimate 

5. Probability of achieving less than 10 percent error in 
production estimate. 

I) . S uppo rtive P rodu ct : o n pat a (for each zone within the region) 

1 . Zone number 

2. Zone name 

3. Zone production estimate 

4. Standard deviation of cone production estimate. 

E. Hist orical Product ion. Data 

1. Number of years of historical data available 

2. Maximum historical production for the region with the 
associated year 
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3. Minimum historical production for the region with the 
associated year 

4.. Mean historical production for the region 

5. Standard deviation of historical production. 

F. Commentary - 160 characters of user input commentary. 

.3.4 Co untry . F.ach LACIF country production estimate report 
shall consist of the following information: 

A. Heading 

1. Title 

2. Country number 

3. Country name 

4 . Crop year 

5. Crop type (spring, winter, or total) 

6. Units indicator 

7. Thresholds applied cn the report data 

• Minimum classification evaluation code 

• Biological stage range 

• Acquisition date range 

• Yield estimate acquisition date range. 

B . Sccu ity Block 

1. CAS monthly report date - user input 

2. Date through which maximum protection must be retained 


3-39 



JSC-10009 


3. " Date through which restricted access must be retained 

4. Date of report generation. 

« 

C. Estimates for Country 

1. Country production estimate 

2. - Standard deviation of production estimate 

3. Coefficient of variation of production estimate 

4. Upper and lower confidence limits of production estimate 

5. Probability of achieving less than 10 percent error in 
production estimate. 

D. Supportive Production Data (for each region within the 
country) 

1. Region number 

2 . Region name 

3. Region production estimate 

4. Standard deviation of region production estimate. 

E. Historical Production Data 

1. Number of years of historical data available 

2. Maximum historical production for the country with the 
associated year 

3. Minimum historical production for the country with the 
associated year 

4. Mean historical production for the country 

5. Standard deviation of historical production. 
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F. Commentary - 160 characters of user input commentary. 

3.5.4 Auxiliary Re p orts . In addition to the area, yield, and 
production estimation reports, auxiliary reports may be optionally 
requested. These auxiliary reports shall consist of the area 
substrata reports, yield strata reports, and WMO position reports. 

3. 5.4.1 Ar ea .Subs trata Reports. The area substrata data shall 
be obtained only for those countries whose base level is sub- 
strata. Each substrata area estimate report shall consist of 
the following information: 

A. Heading 

1. Title 

2. Substrata number 

3. Substrata name 

4. Strata, zone, region, and country to which the substrata 
be longs 

5. Crop year 

6. Crop type (spring or winter or total) 

7. Units indicator 

8. Thresholds applied on the report data 

• Minimum classification evaluation code 

• Biological stage range 

• Acquisition date range. 

1. CAS monthly report date - user input 

2. Date through which maximum protec t ion must be retained 
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3. Date through which restricted access must be retained 

4. Date of report generation. 

C. Area Estimate for Substrata 

1. Substrata area estimate 

'2. Standard deviation of area estimate 

3. Coefficient of variation of area estimate 

4. Upper and lower confidence limits of area estimate 

5. Probability of achieving less than 10 percent error in 
area estimate. 

D. Supportive Area Data 

1. Number of agricultural segments in the substrata 

2. Number of sample segments within the classification 
threshold for aggregation vs. number of sample seg- 
ments allocated to the substrata 

3. Parameters for each sample segment within the substrata 

• Sample segment number 

• Percent classified as wheat 

• Percent other 

• Percent unidentifiable 

• Percent thresholded 

• Classification evaluation code 

• Biological stage 

• Acquisition date(s). 
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E. Historical Area Data 

\ 

1. Number of years of historical data available. 

2 _. Maximum historical area for the substrata with the 

associated year 

3. Minimum historical area for the substrata with the 
associated year 

4. Mean historical area for the substrata 

5. Standard deviation of historical area. 

F. Commentary - 160 characters of user input commentary. 

3. 5.4.2 Yield Strata Reports . Each yield strata estimate report 
shall consist of the following information: 

A. Heading 

1. Date of report generation 

2. Title 

3. Yield strata number 

4. Country to which the yield strata belongs 

5. Crop year 

6. Crop type (spring or winter) 

7. Units indicator. 

B. Yield Strata Data 

1. Yield strata estimate 

2. Variance of yield strata estimate 
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3. Parameters for each area strata within the yield strata 

• Area strata number 

• Area strata name. 

3. 5.4. 3 WMO Posit ion Kcp_oij.s . bach WMO report shall consist of 

t:hc following, information: 

A . Head ing 

1. Date of report generation 

2. Title 

3. Country. 

B. W MO Data (parameters for each WMO station within the country) 

1. WMO station number 

2. Latitude of WMO station 

3. Longitude of WMO station 

4. Yield strata number within which the WMO station is 
located. 

3. 5.4.4 Auxiliary Report Da_t:i Requirements. In or<ier to support 
the auxiliary reports, the following data shall be required in the 
CAS data base: 

A . Y _i e_l_d _ _St r ala Da t a 

• Country ID 

• Yield strata number 

• Crop year 

• Units indicator 

• Yield estimate 
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• Variance of yield estimate 

• Area strata numbers which belong to the yield strata. 
B. : WMO Position Data 


• Country ID 

• WMO station number 

• Latitude 

• Longitude 

• Yield strata number to which the WMO station belongs. 

3.5.5 Report Options . The options for report output are request 
for intermediate aggregation reports and the output device (s) 
for the reports. Reference the display processor (paragraph 2.2.2) 
for further explanation. 

3. 5. 5.1 Intermediate Aggregation Reports . In addition to the 
final aggregation report, all the interim aggregation reports that 
support the final report may be output. 

3. 5. 5. 2 Auxiliary Reports . Each of the auxiliary report types 
may be output for a country per a request at input command time. 
Reference input command processor (paragraph 2.2.1) for further 
explanation. 

3. 5. 5. 3 Output Devices . Any combination of the batch devices or 
CRT may be specified as the output device (s) for report generation. 
The batch devices include a magnetic tape and a line printer. 

3. 5. 5. 4 Data Suppression Flag . At any time during report display, 
the suppression of the output of supportive and historical data 
may be requested or overridden. 
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SECTION 4 


CAS DATA BASE SOFTWARE 


4.1 GENERAL 

Requi remen t s for CAS data base software have been segregated into 
three main functional areas. Paragraph 4.2 addresses requirements 
for data base management services to define, initialize, and update 
CAS data bases in the batch environment. Paragraph 4.3 defines 
requirements for interfacing with CAS applications software (sec- 
tion 3). Paragraph 4.4 defines requirements for interfacing with 
the CAS Data Base Change Program. 

4.2 CAS DATA BASE. MANAGEMENT 


The data base management software shall provide all services nec- 
essary to construct and maintain CAS system data bases. Required 
services shall include capabilities to define, initialize, and 
load CAS data bases and to modify selected data base contents both 
in a batch environment and through, terminal inputs (application) 
processing. All parameters in the data base that contain wheat 
values for area, yield, and production shall he in metric units. 


The following paragraphs address basic, data base 
tions required for data base definition, initial! 
mode processing of raw data inputs received from 
Descriptions of the media, format, and trar.sr. issi 
ternal data source input reflect those defined in 


management fune- 
ral ion, and batch 
external sources, 
on of each ex- 
a upend i x C . 


4.2.1 E xternal Data Sources . An ex t. e rn a 1 
here to .mean any source which, generates da 
at facilities other than the PDP ;i/4f> con 
Building 17 at the JSC. 


a fa source is defined 
inputs for the CAS 
gu rat ion located in 


External data sources are defined as inputs 
lion and Mensuration Subsystem (CAMS), Data 
cessing, and Transmission Subsystem (DAPTS) 
Subsystem (YES), and LACiil allocation data. 


f ror. the C l as s i f i c. a - 
Acquisition, Prepro- 
Yield E s t i na t i on 
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A. CAMS Data , Raw CAMS data inputs shall be received on a 
daily basis for initial processing by a CAS preprocessor 
(see section 5). Subsequent data base updating shall be 
performed by appropriate CAS data base management services 
upon receipt of a preprocessor-generated magnetic tape of 
sample segment classifications ; 

B. DAPTS Data . DAPTS inputs shall be provided to the CAS sys 
tern on a yearly basis. Data provided by DAPTS shall inclu 
spring and winter wheat values for acreage, yield, and pro 
duction by substrata, strata, zone, region, and country as 
applicable to a given country. CAS shall provide storage 
and retrieval facilities for a maximum of (the latest) 15 
years of DAPTS data. DAPTS data shall be provided to CAS 
in' the form of 9-track computer-compatible magnetic tapes. 

C. YES Data . YES -data shall be provided to the CAS on a 
monthly basis. Yield strata-correlated data such as yield 
•strata coefficients and WMO stations as well as meterolog- 

ical data shall be provided to CAS. Up to 15 entries of 
yield estimate data per crop type shall be stored on disk 
for a given crop year. YES data shall be provided to CAS 
in the form of magnetic tape and punched cards. 

D. Allocation Data . Allocation data shall be provided on a 
crop year basis for each LACIE country. Allocation data 
shall define each level which comprises a country hierarch 

4.2.2 Data Base Structuring/Design . Data base structuring/desig 
capabilities shall include provisions for data base security, ini 
tialization, and hierarchical definition of required CAS data 
base(s). The CAS data base shall be sized to hold a minimum of 
2 crop years of data for each data type. 

4. 2. 2.1 Data Base Security . The integrity and security of CAS 
data bases shall be accomplished by operationally controlling 
access to the data base and by utilization of the standard securi 
features offered by the RSX-11D Operating System. Protection of 
programs and data shall be provided by the utilization of a User 
Identification Code (UIC) and passwords. 
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4. 2.2 . 2 Initialization . The process of data base initialization 
includes' those preliminary functions which must occur before a 
data base can be considered available for user access. Functions 
involved in this process may require utilization of certain main- 
tenance capabilities such as data base purging (either a complete 
data base or selective data within a given data base) ; data base 
definition processing; data base loading (e.g., from a previous 
dump or checkpoint tape); and initialization of internal tables, 
indices, pointers, etc. 

All terminal users must be signed off from the affected data base 
prior to initialization processing and, once processing begins, 
shall be locked out from accessing that data base until initializa- 
tion is completed. 

4. 2. 2. 3 Data Base Definition/Hierarchy . The CAS data base shall 
consist of multiple data files. A data file shall be defined for 
each LACIE country for each data type (CAMS, YES, DAPTS) for each 
crop year. Allocation data shall be provided to the CAS each crop 
year for each of the LACIE countries. Allocation data shall define 
the geographical hierarchy employed in CAS data acquisition and 
shall be utilized as the logical structure definition of CAS data 
files. Allocation data shall also be stored on disk for use by 
application and display generation software. 

The data' base management software shall utilize a step-down con- 
catenation of allocation ID’s to perform data storage and retrieval 
operations within a given data file. The actual physical organi- 
zation of a data file shall be transparent to applications inter- 
facing with the data, base management software. 

The following paragraphs address storage and retrieval conventions 
of CAS data files for each data type. 

A. CAMS Data Files . For a given crop year, a data file shall 
be defined for each LACIE country for storage/retrieval of 
CAMS-only data. A physical record within a CAMS data file 
shall be referenced through the concatenation of the follow- 
ing allocation identifiers as applicable to a given country. 

* Country ID 

• Region ID 
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the allocation hierarchy. Each logical record shall con- 
tain values for the following parameters: 

• Epoch year 

9 Crop (harvest) year 

• Spring/winter wheat designator 

• Historical wheat area 

• Historical wheat yield 

• Historical wheat production. 

The storage/retrieval of DAPTS data shall occur at each of 
the following hierarchical levels as applicable to a given 
country. 

• Country 

• Region 

• Zone 

• Strata 

• Substrata. 

C. YES Data Files . A data file shall be defined for each 

LACIE country for storage/retrieval of yield estimate data 
for each crop year. Each physical record of yield estimate 
data shall contain logical records of monthly estimate data 
for a single area stratum, A physical record of yield es- 
timate data shall be referenced through the concatenation 
of the following allocation identifiers as applicable to 
a given country. 

• Country ID 

• Region ID 
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